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Abstract 

A recent analysis of the Fermi Large Area Telescope data provided evidence 
for a high-intensity emission of high-energy gamma rays with a E~ 2 spec- 
trum from two large areas, spanning 50° above and below the Galactic cen- 
tre (the "Fermi bubbles"). A hadronic mechanism was proposed for this 
gamma-ray emission making the Fermi bubbles promising source candidates 
of high-energy neutrino emission. In this work Monte Carlo simulations re- 
garding the detectability of high-energy neutrinos from the Fermi bubbles 
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with the future multi-km 3 neutrino telescope KM3NeT in the Mediterranean 
Sea are presented. Under the hypothesis that the gamma-ray emission is com- 
pletely due to hadronic processes, the results indicate that neutrinos from 
the bubbles could be discovered in about one year of operation, for a neu- 
trino spectrum with a cutoff at 100 TeV and a detector with about 6 km 3 of 
instrumented volume. The effect of a possible lower cutoff is also considered. 
Keywords: neutrino telescope, Fermi Bubbles, KM3NeT 

1. Introduction 

In the last decade a new generation of telescopes revealed a large variety 
of astrophysical high-energy gamma-ray sources. Telescopes such as Fermi 
[U, HESS 0, VERITAS [3] and MAGIC [4] can identify sources of gamma 
rays with energies from about 20 MeV to about 100 TeV with sufficient angu- 
lar resolution to study the source morphology of extended Galactic sources. 
They can also measure the gamma-ray energy spectrum with high precision. 
From these measurements the presence of cosmic acceleration processes has 
been confirmed in a large variety of known and unknown galactic and extra- 
galactic sources [5j [6j . 

For the full understanding of the underlying mechanisms other probes 
that can complement the information from gamma-ray detection are needed. 
As neutrinos are mainly produced via proton-proton (jyp) and proton-gamma 
(p7) interactions, they are an unambiguous signature of hadronic accelera- 
tion. 

Neutrino telescopes have been constructed, and more sensitive ones are 
proposed, to detect these high-energy neutrinos. Neutrino flux estimates 
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indicate that detectors of km 3 -scale instrumented volume are required [7J. 
The IceCube [8] detector located at the South Pole, about 1 km 3 in size, is 
taking data in its final configuration since 2011. In the northern hemisphere 
ANTARES is the largest operating detector located in the Mediterranean Sea 
40 km off the French coast near Toulon [9]. It has an instrumented volume 
of about 0.01km 3 and it is taking data in its final configuration since 2008. 
The KM3NeT consortium [10J proposes the construction of a multi-km 3 tele- 
scope in the Mediterranean Sea (KM3NeT). From their respective locations, 
KM3NeT and ANTARES will detect upward-going neutrinos from about 
3.57rsr of the sky, including the Galactic centre and most of the Galactic 
plane, where many TeV gamma-ray sources are located [SJ [7J . 

In this work the response of the KM3NeT telescope to neutrinos from the 
Fermi bubbles (see below) is studied with Monte Carlo simulations under 
the hypothesis that a hadronic mechanism is responsible for the gamma- 
ray emission. In particular, the neutrino flux that could be discovered by 
KM3NeT is determined as a function of the number of observation years. 

2. The Fermi bubbles 

A recent analysis of Fermi-LAT data [TTJ revealed an intense gamma-ray 
emission from two large areas above and below the Galactic centre. The 
detected gamma-ray emission has the following characteristics: 

• The emission areas are symmetric with respect to the Galactic plane 
and extend up to 50 degrees (lOkpc) north and south the Galactic 
centre, with a width of 40 degrees in Galactic longitude. 
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• The gamma energy spectrum, measured from ~ 1 GeV to ~ 100 GeV, 
is compatible with a power-law spectrum described by E 2 d^^/dE « 
3-6 x 10- 7 GeVcm- 2 s- 1 sr- 1 . 

• The emission is homogeneous within the bubbles. No significant differ- 
ences were found between the northern and the southern bubble. 

• The edges of the bubbles seem to be correlated with ROSAT X-ray 
maps at 1.5— 2keV, while the inner parts are correlated with the hard- 
spectrum microwave excess known as WMAP haze [HJ H2J [13] . 

Several mechanisms have been suggested in the literature to explain 
this gamma-ray emission. Most of these explanations rely on leptonic pro- 
cesses and include inverse Compton scattering by electrons, either produced 
by long-lasting energy injections near the Galactic centre [11 J or acceler- 
ated through the second-order Fermi process generated by turbulent plasma 
throughout the entire bubble |14j . Other processes such as recent transient 
AGN activity near the Galactic centre [15] , millisecond pulsars [16] and dark 
matter annihilation [17] have also been put forward. 

Currently the observed features cannot be fully explained by these pre- 
dominantly leptonic processes, which has led to a proposal for an underlying 
hadronic process [18]. A cosmic ray population associated with long-time- 
scale star formation in the Galactic centre (of the order of 10 10 years) was 
injected into the bubbles where it interacts with the ambient matter and 
produces high-energy gamma rays through 7r° decay. Pair production is re- 
sponsible for the generation of secondary electrons producing the synchrotron 
radiation observed as microwaves. In this scenario neutrinos are produced 
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through the decay of charged pions. The leptonic and hadronic mechanisms 
require very different galactic diffusion characteristics. 

Recently the morphological and spectral characteristics of the WMAP 
haze have been reviewed [T3], with important implications for the possible 
production scenarios. The author considers all the proposed models to ex- 
plain the gamma and the WMAP haze emission and concludes that hybrid 
emission scenarios will likely be required. 

The detection of high-energy neutrinos can distinguish between hadronic 
and leptonic models. Only if a hadronic process is totally or partially respon- 
sible for the production of gamma rays from the Fermi bubbles, neutrinos 
will also be produced. The KM3NeT telescope, to be located in the deep 
Mediterranean Sea, will be the ideal instrument for the observation of neu- 
trinos from the Fermi bubbles. 

3. The detector 

The KM3NeT consortium aims at installing a deep-sea research infras- 
tructure hosting a multi-cubic-kilometre high-energy neutrino detector in the 
Mediterranean Sea. During the Design StudjQ several technical solutions 
were investigated and the results reported in the Technical Design Report 
|20j . During a Preparatory Phase^J a final design concept was defined. Cur- 
rently prototyping activities are in progress to prepare for the construction 
phase. 

The detection principle relies on the detection of Cherenkov light induced 

Supported by the EU in FP6, contract n° 011937 

2 Supported by the EU in FP7, grant agreement n° 212525 
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by secondary charged particles produced in neutrino-nucleus interactions in- 
side or near the detector. The Cherenkov light is detected by photomultiplier 
tubes (PMTs) contained in glass spheres that are designed to resist the hy- 
drostatic pressure of the deep-sea environment. These instrumented spheres 
are called Optical Modules (OMs). 

The geometry of the detector [21] simulated in this work consists of a 
three-dimensional array of OMs attached to vertical structures (Detection 
Units, DUs). An array of DUs constitutes a detector building block. Several 
building blocks form the full detector of about 300 DUs. A DU, which is 
anchored to the sea floor and kept upright by a submerged buoy, consists of 
horizontal bars equipped with 2 OMs, one at each end. Adjacent bars are 
oriented orthogonally to each other. The DUs are connected to shore by an 
electro-optical cable. Each OM consists of 31 3-inch PMTs [22] housed inside 
a 17-inch pressure-resistant glass sphere covering almost a 4tt field of view. 

For this work, a detector of 308 DUs (two building blocks of 154 DUs each, 
arranged uniformly in a circular area) with an average separation between 
adjacent DUs of 180 m has been simulated. Each simulated DU consists of 
20 bars of 6 m length with a vertical spacing of 40m (40 OMs per DU). The 
total instrumented detector volume is about 6 km 3 . 

4. Monte Carlo simulation 

In order to study the sensitivity of the KM3NeT telescope for the de- 
tection of neutrinos from the Fermi bubbles, we use a Monte Carlo simula- 
tion framework. The framework is based on the ANTARES software [23J, 
modified for a km 3 -scale detector and using OM properties appropriate for 
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KM3NeT. The simulation chain consists of the generation of muon neutrinos 
from the bubbles, the generation of atmospheric muon and neutrino back- 
grounds, the neutrino charged-current interactions, the propagation of the 
produced muons in rock and sea water, the generation of Cherenkov light, 
the 40 K background and the digitisation of the PMT signals. Optical prop- 
erties of the sea water and the PMT characteristics are taken into account in 
the simulation. The depth and the optical water properties measured at the 
Sicilian Capo Passero site have been used [20] • Background light due to the 
presence of 40 K in salt water has been simulated adding an uncorrelated hit 
rate of 5 kHz per PMT and a time-correlated hit rate of 500 Hz per OM (two 
coincident hits in different PMTs inside the same OM). These parameters 
have been estimated with a complete simulation based on GEANT4 |24j. 
The 5 kHz of uncorrelated hit rate is consistent with the baseline of 50 kHz 
measured with the 10-inch PMTs of the ANTARES experiment [25J. 




The muon track direction is reconstructed from the simulated arrival 
times of Cherenkov photons and the PMT positions. The reconstruction 
algorithm is based on a hit selection and on a maximum likelihood fit that 
uses probability density functions for the photon arrival times at the PMTs. 
In addition to the positions and track directions, the number of hits used for 
the reconstruction (iVhit) and a track fit quality parameter (A) are given as 
output. The A parameter is determined from the likelihood and from the 
number of compatible track solutions found by the algorithm and is used to 
reject badly reconstructed events [26]. The A^it parameter is correlated with 
the muon energy. 
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4-1. Neutrinos from the Fermi bubbles 

Muon neutrinos from the Fermi bubbles were generated homogeneously 
within two circular regions of 19° radius around two positions in the sky at 
the equatorial coordinates declination 5 = — 15° and right ascension a = 243° 
for the northern bubble and 8 = —44° and a = 298° for the southern bubble. 
The simulated neutrino energy is between 10 2 and 10 s GeV. 

Figure [T] shows the percentage of Monte Carlo neutrinos that are up-going 
in local detector coordinates. Since only up-going events can be unambigu- 
ously classified as neutrino candidates, Fig. [T] represents the Fermi bubbles' 
visibility for a neutrino telescope in the Mediterranean Sea. Black points are 
the edges of the two bubbles as reported in [11] . The average visibility for a 
detector located at a latitude of 36° 16' N is 58% of the time for the northern 
bubble and 80% for the southern bubble. We note that for IceCube only a 
small fraction of the Fermi bubbles lies below the horizon (solid red line in 
Fig. [I]), leading to a significantly reduced sensitivity compared to a similar 
detector in the Mediterranean Sea [2*7] . 

The generated Monte Carlo events can be weighted to reproduce different 
assumed neutrino spectra. The general energy dependence of the neutrino 
flux used in this analysis is a power law spectrum with a spectral index of 
—2 and an exponential cutoff: 

d ^t = K Qv .E-*.e- E ' E % (1) 

where E c is the cutoff energy. For this analysis we consider three cases: 
d&^/dE, i.e. a pure power law spectrum with spectral index —2, as well 
as d$l 00 /dE and d&f/dE with cutoffs at 100 TeV and 30 TeV, respectively. 
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The cutoff at E c = 100 TeV is consistent with an assumed cutoff of the proton 
energy distribution in the Fermi bubbles at a few PeV (corresponding to the 
position of the "knee" in the cosmic ray energy spectrum). The second cutoff 
value considered, E c = 30 TeV, is more conservative and does not require the 
assumption that the Fermi bubbles populate the cosmic ray spectrum up to 
the knee. 




Galactic longitude (degrees) 

Figure 1 : Location of the Fermi bubbles as simulated in the analysis in Galactic coordi- 
nates. The color code indicates the fraction of Monte Carlo neutrinos that are up-going 
in local detector coordinates. The measured bubble edges ([TT] Table 1) are also reported 
(black points). The red line represents the separation line between the northern and 
southern hemisphere. 

Under the hypothesis that the source is transparent to gamma rays and 
that the mechanism responsible for the gamma-ray emission is hadronic, the 
neutrino spectrum was estimated from the measured gamma-ray spectrum 
following the prescription described in [19]. The measured Fermi bubble 
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gamma-ray spectrum is consistent with a E~ 2 shape and a normalisation 
i^Q 7 « 4 x 10~ 7 GeV cm -2 s _1 sr^ 1 . Using this gamma-ray flux and taking 
into account the 0.7 sr solid angle of the two simulated bubbles, we estimated 
the corresponding muon neutrino plus antineutrino flux to have Kq v ~ 1 x 
1(T 7 GeVcm^s- 1 . 

4-2. Atmospheric muons 

Cosmic rays entering the atmosphere produce extensive air showers that 
contain high-energy muons. Although the sea water above the detector serves 
as a shield, many of these muons reach the detector. Therefore down-going 
tracks are excluded from the analysis. However, atmospheric muons mis- 
reconstructed as up-going remain a significant background. Simulating at- 
mospheric muon background in a large detector requires huge amounts of 
CPU time and data storage. Atmospheric muons were generated with the 
fast MUPAGE code This code, which is based on a full Monte Carlo 
simulation of primary cosmic ray interactions and shower propagation in the 
atmosphere, provides a parameterisation of the underwater flux of atmo- 
spheric muons including also multi-muon events ("muon bundles"). 

Atmospheric muons were generated in the range 1 TeV < Eb < lOTeV, 
where Ej, is the sum of the energies of all single muons in the bundle. This 
sample is statistically equivalent to a live time of 2 hours. To efficiently 
increase statistics in the high-energy region an additional sample equivalent 
to 8 days of live time was generated with Eb > 10 TeV. The Monte Carlo 
samples were then reweighted to the relevant experimental live time. 
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4-3. Atmospheric neutrinos 

A large number of secondary particles is produced in a cosmic ray in- 
teraction in the atmosphere. The production of pions and kaons and their 
subsequent decay chains produce a large flux of atmospheric neutrinos, called 
"conventional flux" , with an energy distribution proportional to -E~ 3 ' 7 at neu- 
trino energies E v larger than about 1 TeV. The production of neutrinos from 
heavy quark is responsible for a high-energy component [E v > 10 TeV) in the 
atmospheric neutrino spectrum. This contribution, called "prompt", is not 
well known, but several predictions are available [29, 30J. The atmospheric 
neutrino flux is an irreducible background for the detection of neutrinos of 
cosmic origin. 

The atmospheric muon neutrino and antineutrino background was gen- 
erated in the energy range 10 2 GeV < E v < 10 8 GeV and over the full solid 
angle. The events are weighted to reproduce the conventional atmospheric 
neutrino flux following the Bartol model [31] . The uncertainty on the Bartol 
normalisation factor is taken to be about 25-30% [32] ■ A prompt contribu- 
tion is also taken into account through the event weights. The models in 
[29l [30] have been considered and the model with the highest neutrino flux, 
corresponding to the highest prediction of the Recombination Quark Parton 
Model (RQPM) |33J, has been used in the present analysis. 

5. Flux discovery potential 

The present analysis was restricted to events reconstructed as up-going 
and located within 19° around the centre of each Fermi bubble. The bubbles 
are considered as discovered if the number of detected events in a given 
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detector live time has a probability of a = 2.85 x or less to originate 
purely from background in 1 — (3 = 50% of all experiments. This corresponds 
to a significance of 5a (area of the one-sided Gaussian tail). 

The signal flux required to claim a discovery is calculated from the sim- 
ulated average number of background events, (nb ac k)- First, the minimum 
(critical) number of events, no, that satisfies 

oo 

22 ^(^obsK^back)) < a (2) 

™obs = ™0 

is determined, where -P(n b s |(nb a ck)) is the Poisson probability for observing 
?2 bs events given (nb ac k)- The value no is the minimum number of events 
required to claim a deviation from the background-only hypothesis with a 
statistical significance defined by the p-value a. The confidence level CL. is 
related to p via CL. = 1 — p. 

In case of the presence of a signal of strength n a , the probability to observe 
n or more events is related to the statistical power (1 — (3) by 

oo 

^ P(n ohs \(n a ) + (™back)) = 1-/3. (3) 
™obs=™0 

The strength of such a signal would lead to an observation with a p-value 
less than a in a fraction {1 — 0) of the experiments. 

In order to determine the minimum flux <p a needed for a discovery (dis- 
covery flux), the cuts on A and iVhit were varied and the minimum in the 
Model Discovery Potential (MDP) [S], defined as MDP = n Q ((n back ))/(n s ) 
was sought. Here, (n s ) is the number of signal events resulting from the 
model flux (fij) with K Gv w 1 x 10~ 7 GeV cm -2 s _1 after cuts. The discovery 
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flux is then related to the MDP by 
Q = 0o • MDP 



n a { (n hack )) 

(n s ) 



(4) 



Figure [2] shows the cumulative distributions of the numbers of events per 
year, reconstructed as up-going in the bubble region, for the different Monte 
Carlo samples as a function of A and Nat- Figure [2] demonstrates that most of 
the events due to atmospheric muons that are misreconstructed as up-going 
can be rejected by imposing an appropriate cut on A, while an appropriate 
cut on A^hit helps in the rejection of atmospheric neutrino background events 
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Figure 2: Cumulative distribution of the number of events per year reconstructed as 
up-going in the bubble region as a function of A (left panel) and of iVhit (right panel) for 
the different Monte Carlo samples. For the right panel a cut A > — 10 was applied. 

The cuts required to minimise the MDP vary for the different flux assump- 
tion and are given in Tab. [1} Figure [3] shows the spectrum of reconstructed 
signal and background up-going events per year in the region of the bubbles 
as a function of the simulated neutrino energy. The signal energy spectrum 
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Figure 3: Number of events per year reconstructed as up-going in the bubble region as 
a function of the simulated neutrino energy, for the signal and the atmospheric neutrino 
background with and without the cuts required to minimise the MDP (5a C.L., 50% 
probability). The assumed normalisation factor is Kq v = 1 x 10~ 7 GeVcm~ 2 s _1 . The bin 
width in log 10 (E v ) is 0.1. 

from the bubbles, after the cuts, is centred at about 500 TeV for / dE 
(see Eq. Q) (a), at about 40 TeV for d$l°°/dE (b) and at about 25 TeV for 
d&l°/dE (c). The numbers of events at the reconstruction level and after 
the cuts in one year of data taking are reported in Tab. [T]for a source with 
a flux normalisation factor of K 0u = 1 x 10 -7 GeVcm~ 2 s _1 . 

With the small simulated live time of the atmospheric muon samples the 
number of events per year at cut level due to the atmospheric muons is zero 
(see Tab. [T]). By extrapolating in these data sets the trends of the numbers 
of events per year as a function of A and iV hit (Fig. [2]), it was estimated that, 
in the cases of neutrino spectra with cutoffs, the contribution of atmospheric 
muons to the background does not exceed 25% of the atmospheric neutrino 
background events reported in Tab. [T} A 25% increase in the number of 
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Table 1: Expected numbers of events in one year of data taking for the three samples 
of simulated events: neutrinos from the Fermi bubbles (z4ig)j neutrinos from atmospheric 
background (conventional plus prompt RQPM) (V a tm) and muons from atmospheric back- 
ground (na,tm)- The cut values (A and iVhit) and the numbers of events at the reconstruc- 
tion level (reco level) and after the cuts (cut level) also also given. The cuts are those 
optimised for a discovery at 5ct C.L., 50% probability, for a source with a normalisation 
factor of Kqu = 1 X l(T 7 GcVcm- 2 s _1 . 

background events corresponds to an increase of 10% in the flux limits. In 
the case of a pure power-law spectrum the contribution of atmospheric muons 
to the background is negligible. 

In Fig. [i] the discovery fluxes at 5a C.L., 50% probability and 3cr C.L., 
50% probability for the neutrino spectra considered (see Eq. Q) are shown 
as functions of the observation time. The variation in the discovery fluxes 
due to the uncertainties in the atmospheric neutrino flux have been evaluated 
taking into account only the uncertainty on the normalisation factor of the 
conventional Bartol flux. The RQPM prompt model component contributes 
with about 46% and 10% in the number of neutrino background events after 
the cuts quoted in Tab. fl] for the pure E~ 2 spectrum and for the exponential 
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cutoff spectra, respectively. A 25% uncertainty in the conventional Bartol 
atmospheric neutrino flux corresponds to a variation in the discovery flux of 
about 7% for d$™/dE (see Eq. Q) and 10% for d® 1 ™ / dE and d&f/dE. 
These variations are shown in Fig. [4] as bands around the discovery-flux 
curves. 

If the neutrino normalisation factor is of the order of 1 x 10~ 7 GeV cm -2 s _1 
the discovery of neutrinos from the Fermi bubbles is expected in about one 
year of data taking for a E~ 2 neutrino spectrum with a cutoff at 100 TeV. 
The first evidence (3a C.L., 50% probability) could be obtained after a few 
months of data taking. For the more severe cutoff at 30 TeV the discovery 
is predicted to be achieved in about 5.5 years and first evidence in about 
2.5 years. 

Recently an analysis has been performed on 43 months of Fermi data in 
the energy range 20 GeV < E 1 < 300 GeV [36J . This analysis covered several 
regions near the Galactic centre, some of which overlap with the bubble re- 
gions. This revealed the presence of an excess in the gamma ray flux around 
130 GeV. This peak has been interpreted as an indication of dark matter 
annihilation into two photons. An alternative explanation, put forward in a 
separate analysis of the same data (in the range 80 GeV < E 1 < 200 GeV) 
[37] . is that the excess can be interpreted as a steepening of the power law 
spectrum of the Fermi bubble gamma rays at around 130 GeV. Such a steep- 
ening has a significant effect on the results of the present analysis. There- 
fore, a Fermi bubble neutrino spectrum that follows a power law spectrum 
according to d<k™/dE with K 0u (E) = 10~ 7 GeV cm" 2 s" 1 up to 130 GeV and 
extends to higher energies proportional to E~ 2 ' 3 , has been investigated. In 
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Figure 4: Discovery fluxes as functions of the observation years for 5er C.L., 50% prob- 
ability and 3fi C.L., 50% probability, for the three neutrino spectra assumed. The bands 
represent the variation due to the uncertainty on the normalisation factor of the conven- 
tional Bartol neutrino flux. The long-dashed line indicates the estimated neutrino flux 
value. 

this case the resulting live time required for discovery is 6 years, which is 
well in the expected lifetime of the KM3NeT detector. 

Since neutrino telescopes have a modular design, the detector concept 
allows for a staged implementation with continuously increasing science ca- 
pabilities. Therefore, the capability of the telescope to detect neutrinos from 
the Fermi bubbles has been explored as a function of an increasing number 
of DUs. The discovery fluxes for telescope stages composed of 50, 100 and 
154 DUs have been calculated and are shown in Fig. [5] for 5cr C.L., 50% 
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Figure 5: Discovery fluxes as a function of the observation years for a 5a C.L., 50% 
probability (a) and 3tr C.L., 50% probability (b) for a E~ 2 neutrino spectrum with a 
100 TeV cutoff for different detector sizes (see text). The bands represent the variation 
due to the uncertainty on the normalisation factor of the conventional Bartol neutrino 
flux. The long-dashed lines indicate the estimated neutrino flux value. 

probability (a) and 3a C.L., 50% probability (b). The detector geometries 
have the same average DU distance (180 m) and provide increasing active 
volumes of 1, 2.1 and 3.2 km 3 for 50, 100 and 154 DUs, respectively. The 
number of years needed for the discovery does not decrease linearly with the 
increasing number of DUs, being about 8 years for 50 DUs, 4 years for 100 
DUs and 2.5 years for 154 DUs. A neutrino spectrum with normalisation 
Kq u ~ 1 x 10~ 7 GeV cm" 2 s -1 with a 100 TeV cutoff energy is assumed here. 

6. Conclusions 

The discovery of an intense gamma-ray flux from two large areas around 
the Galactic centre, the Fermi bubbles, has stimulated estimates of the ca- 
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pability of neutrino telescopes to discover neutrinos from this region [38, 39J. 
The Fermi bubbles extend over a large part of the sky and have a significant 
intensity in gamma rays. Assuming these gamma rays are of hadronic ori- 
gin and the spectrum extends to the multi-TeV range, this analysis shows 
that the bubbles will be the most intense high-energy neutrino source near 
the Galactic centre. Telescopes located in the Mediterranean Sea, such as 
ANTARES and KM3NeT, have a large visibility of the region around the 
Galactic centre. Due to its multi-km 3 scale, the future KM3NeT telescope 
is the ideal instrument to observe neutrinos from the Fermi bubbles. Note 
that, due to its location, IceCube has a considerably reduced visibility of the 
bubbles. 

In this paper we present an analysis addressing the sensitivity of KM3NeT 
to a neutrino flux from the Fermi bubbles that is consistent with the mea- 
sured gamma ray flux, assuming that the latter fully originates from hadronic 
processes. Based on a complete chain of Monte Carlo simulations, the time 
of KM3NeT operation that is required to detect this neutrino flux at a sig- 
nificance of 5a with 50% probability has been estimated. The result depends 
strongly on the shape of the neutrino energy spectrum assumed. In the case 
of an E~ 2 spectrum with an exponential cutoff at 100 TeV, we expect to 
obtain a 3cx-evidence in a few months and to claim a discovery after about 
1 year of data taking. The non-observation of a signal would severely con- 
strain models of neutrino production through hadronic acceleration processes 
in the Fermi bubbles. 
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